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This article reports a measurement of the cross section of prompt isolated photon pair production in pp collisions 
at a total energy y/s = 1.96 TeV using data of 5.36 fb~^ integrated luminosity collected with the CDF II detector 
at the Fermilab Tevatron. The measured cross section, differential in basic kinematic variables, is compared with 
three perturbative QCD predictions, a leading order (LO) parton shower calculation and two next-to-leading 
order (NLO) calculations. The NLO calculations reproduce most aspects of the data. By including photon 
radiation from quarks before and after hard scattering, the parton shower prediction becomes competitive with 
the NLO predictions. 
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C/^ , 1. Introduction 
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. The measurement of the production cross section of two energetic isolated central photons (diphotons) in high 

' energy hadron collisions is important for testing standard model (SM) predictions in the domain of searches 

^ for undiscovered particles and new physics. Photons originating from hard collisions of hadrons ("direct" 
QJ or "prompt" photons) are an ideal probe for testing perturbative Quantum ChromoDynamics (pQCD) and 
JL , soft-giuon resummation methods implemented in theoretical calculations because they do not interact with 
other final state particles, and their energies and directions can be measured with high precision in modern 
] electromagnetic calorimeters. Prompt diphoton production with large invariant mass creates an irreducible 
'— 'i background in searches for a low mass Higgs boson decaying into a photon pair as well as in searches for 
^ ' new phenomena, such as new heavy resonances |^ , extra spatial dimensions [1, 0] or cascade decays of heavy new 
^ , particles [Ij. Precise measurements of the diphoton production differential cross sections for various kinematic 
J — ■ variables and their theoretical understanding are thus very important for these searches. 

C^l . The basic mechanisms of prompt diphoton production in hadron collisions are quark- antiquark annihilation 
OO ' qij — ^ 77 quark-gluon scattering gq — ^ 77(7, and gluon-gluon fusion gg 77. The respective basic diagrams 
arc shown in Fig. [TJ At the Tevatron, the dominant mechanism is quark-antiquark annihilation. In quark- 
^\ gluon scattering, most of the time at least one of the two photons is emitted almost parallel to the scattered 
\ quark. Contributions from this mechanism are therefore suppressed by requiring isolated prompt photons. 
^""^ . Each mechanism can be modeled by calculating the respective matrix element for the specific event kinematics. 
Matrix element calculations of Leading Order (LO) in the strong coupling are relatively simple and are thus 
implemented in advanced parton shower Monte Carlo (MC) event generators which allow for gluon and 

photon radiation as well as multiple interactions in the colliding beams. By including radiation before and after 
^ the hard scattering, parton shower generators take into account soft gluon and photon emissions, thus resulting 
in an effective resummation of all of the Leading Logarithmic (LL) terms in the cross section to all orders of the 
strong and electromagnetic couplings constants. Next-to-Leading Order (NLO) calculations [9l-[Tl| additionally 
include one-loop corrections at the cost of not featuring realistic multi-particle event representations as the LO 
generators do. Recent NLO calculations include an analytical resummation of the cross section for initial-state 
gluon radiation to all orders in the strong coupling constant [llj, reaching a higher logarithmic accuracy than 
in the parton shower Monte Carlo generators. By this method, all soft gluon emissions in the initial state are 
taken into account, and reliable predictions for the low diphoton transverse momentum region are possible. A 
fixed-order NLO calculation implemented by the diphox program also accounts for the case where a final 
state quark loses almost all of its energy to the photon detected in the event [diagram (e) of Fig. [Tj [12]. This 
process is called "fragmentation" and, in contrast to final state photon radiation in parton showering, it involves 
non-perturbative calculations. One or both photons in the event may come from fragmentation. The case where 
both photons come from fragmentation of a single quark is also possible, but is not included in calculations, as 
in this case the photons are nearly collinear and thus non-isolated most of the time. 

Diphoton measurements have been previously conducted at fixed-target [Tsj and collider experiments [T3 - 
[iTj . The most recent measurements [l^ [13 were compared with the same pQCD calculations examined in 
the present work and large discrepancies were found between the data and a LO matrix element calculation 
supplemented with a parton shower model, suitable for simulation of the backgrounds in searches of a low mass 
Higgs boson and of new phenomena. This work shows that the inclusion of photons radiated from initial and 
final state quarks drastically improves the comparison of the parton shower calculation with the data. 
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Figure 1: Basic diagrams for prompt diphoton production: (a-b) direct, (c-d) one-photon radiation from an initial- (ISR) 
or final-state quark (FSR), (e) fragmentation where one photon is emitted along the direction of a final-state quark 
taking almost all of its energy. The symbol ® denotes the non-perturbative mechanism of the fragmentation process 
(FRAG). 



2. Measurement 

The reported measurement was conducted using data of total integrated luminosity 5.36 fb~^ collected with 
the Collider Detector at Fermilab (CDF) [3l at the Tevatron "pp collider. CDF is composed of a central 
spectrometer inside a 1.4 T magnetic field, surrounded by electromagnetic and hadronic calorimeters and muon 
detection chambers. The inner spectrometer measures charged particle tracks with a transverse momentum 
(pt) precision of Apx/Px — 0.07%(GeV/c)^^. The central calorimeters cover the region \t]\ < 1.1, with 
an electromagnetic (hadronic) energy resolution of (t{Et^)/Ej' — 13.5%/-\/^r(GeV) ® 1.5% {a{E'T)/ET = 
50%/-\/&r(GeV) ® 3%) and a tower segmentation of A77 x A(j) ~ 0.1 x 15°. Photons are reconstructed in 
clusters of up to three towers [l^. criteria are imposed on the profile of the shower to match expected 
patterns. Two main cuts are applied; (i) the photon transverse energy is required to be Et > 17 GeV for the 
first photon in the event and Et > 15 GeV for the second photon; (ii) the calorimeter isolation energy in the 
isolation cone around each photon is required to be less than 2 GeV. 

The background from 7-l-jet and dijet events, where one or two jets are faking a photon, is subtracted 
with a method using the track isolation as the discriminant between signal and background [2l|. It is based 
on the substantial difference of the track isolation distribution for signal photons (nearly exponential) and 
for background photons (nearly flat). The advantages of this method are that (i) it has little sensitivity to 
multiple interactions in the colliding beams, so that the signal-background separation does not degrade at high 
instantaneous luminosity, and (ii) it has high efficiency and good track momentum resolution, implying minimal 
degradation of the signal-background separation due to instrumental effects. The signal fraction is determined 
by summing the probabilities of an event to be pure signal, pure background, or mixed photon pair. These 
probabilities are obtained by solving a 4 x 4 matrix equation using the observation value (0 or 1) for all four 
combinations of the leading or sub-leading photon having track isolation below or above 1 GeV/c as an input. 
The matrix is constructed from the i?T-dependent efficiencies of signal and background photons passing the 
track isolation cut. A threshold cut of 1 GeV/c is determined by maximizing the separation between signal 
and background. The efficiencies are determined from Monte Carlo (MC) 7+jet and dijet samples, which are 
produced using the PYTHIA event generator Q. PYTHIA events are fully simulated thr oug h the detector and 
trigger and are reconstructed with the CDF II simulation and reconstruction software [22|. With this matrix 
technique the full correlations between the two photons in the event are properly taken into account. Tests are 
made for underlying event contributions in complementary cones to the photon reconstruction cone [2^ 1 and 
also using isolated tracks in dijet events. The systematic uncertainty in the signal fraction with this method is 
of the order of 15-20%. 

The diphoton production cross section differential in a kinematic variable is obtained from the histogram of 
the estimated signal in the selected variable. The average cross section in a bin of the variable is determined by 
dividing the bin content by the trigger efficiency, the diphoton selection efficiency and acceptance, the integrated 
luminosity and the bin size. The diphoton trigger efficiency is derived from data [ij. It is consistent with 100% 
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over all of the kinematic range with a flat uncertainty of 3%. The selection efficiency is determined from data 
and MC with an iterative method. In the first pass the efficiency is determined from a fully simulated and 
reconstructed pythia diphoton MC sample by dividing the number of events passing all selection cuts by the 
number of events passing only the kinematic cuts on the photon £^t, angular separation and isolation at 
the event generation level. The efficiency denominator is corrected for the "underlying event" from collision 
remnants which make the efficiency obtained from pythia too high by removing events from the denominator 
through the isolation cut. This correction is derived by running pythia with and without underlying event 
and amounts to a constant factor of 0.88 per event. A flat 6% uncertainty in the selection efficiency (3% per 
photon) accounts for possible inaccuracies in the pythia model for the underlying event. The signal events of 
the data are corrected for the preliminary efficiency. The data are then used to reweight the pythia events 
and obtain a more accurate representation of the true diphoton distribution. The efficiency is determined using 
the reweighted pythia sample and corrected for luminosity dependence, derived from a comparison of the 
vertex multiplicity distribution in data and pythia MC — e+e" events. The systematic uncertainty in the 
efficiency resulting from the luminosity dependent correction grows linearly from 1.8% for < 40 GeV to 3% 
at Et = 80 GeV and remains constant above this point. Finally, a 6% constant uncertainty comes from the 
Tevatron integrated luminosity (2^ . 

The — T' e+e~ sample is used for calibration by applying a "diphoton-like" event selection, i.e. by imposing a 
diphoton selection with the same trigger but allowing for a track associated with each of the two electromagnetic 
objects in the event. The electrornagnetic energy scale in data and MC is corrected by tuning the Z^ — )■ e+e~ 
mass peak to the world average [25| and a systematic uncertainty from this correction is estimated to grow 
linearly from at Et < 40 GeV up to 1.5% at E^ — 80 GeV and remain constant above this point. The 
difference in the photon identification efficiency between data and MC is estimated from the Z° e+e^ sample 
[l| and added as a systematic uncertainty to the measurement. All systematic uncertainties in the cross section 
measurement are added in quadrature. 



3. Predictions 

The results of this measurement are compared with three theoretical calculations: (i) the fixed NLO predic- 
tions of the DIPHOX program [^] including parton fragmentations into photons 12], (ii) the predictions of the 
RESBOS program [ll| where the cross section is accurate to NLO, but also has an analytical initial state soft gluon 
resummation, and (iii) the predictions of the pythia program 6] which features a realistic representation of the 
physics events by including parton showering. Initial (ISR) and Final State Radiation (FSR) and an underlying 
event model. Diphoton events were selected from an inclusive 7-I-X pythia sample (X=7 or jet), thus including 
the qq — >■ 77 and gg 77 processes (56%) as well as the qq 577isn i 91 ^TTisn ^-'^d gq — qn-ps-a processes 
(44%). This type of calculation effectively resums the cross section for gluon and photon radiation both in the 
initial and the final state. Fig. [5] shows the individual contributions to the cross section as a function of the 
diphoton invariant mass, transverse momentum and azimuthal difference. Initial-state radiation (ISR) photons, 
in particular, produce substantially different distributions than ME and final-state radiation (FSR) photons, 
having a harder transverse momentum spectrum and stronger low-A0 tail in the azimuthal difference spectrum. 
In leading order, this can be attributed to the fact that FSR occurs in quark-gluon scattering [diagram (d) of 
Fig. [1], whereas ISR occurs both in qq annihilation [diagram (c) of Fig. [1] and quark-gluon scattering, and 
the luminosity of quark-gluon states falls off more rapidly with the parton momenta than the luminosity of 
qq states [lH. Double radiation processes in minimum bias dijet events, such as qq — )■ 'Z'Z7isr/fsr7isr/fsr' 
qq gg7isR/FSR7isR/FSR' 59 ^ ff97isR/FSR7isR/FSR' m ff57isR7isR and gg <7g7FSR7FSR' were also examined 
but their overall contribution was estimated to only ~ 3% of the total, having no significant effect to any 
kinematical distribution. Therefore, these processes were not included in the pythia calculation. 

All calculations are subject to the experimental kinematic and isolation cuts, diphox accounts for the 
gg — >■ 77 process in LO only. The predictions of RESBOS are restricted to the invariant mass range from 2mb — 9 
GeV/c^ to 2mt — 350 GeV/c^, where rrib and rut are the masses of the bottom and top quarks, respectively. 
NLO theoretical uncertainties are estimated by varying the fragmentation (in diphox only), renormalization 
and factorization scales up and down by a factor of two relative to the default scale /x = M/2 of diphox 
and /i = M of RESBOS, and for the NLO PDF uncertainties (in both DIPHOX and RESBOS) by using the 20 
CTEQ6.1M eigenvectors [26*1. Theoretical uncertainties are shown only in the plots of the relative deviations 
of the data from each calculation, in the form (data— theory) /theory, as functions of the selected variable. No 
relative deviations are shown for the pythia 77 calculation. The benchmark parton showering MC calculation, 
compared in detail with the data, is pythia 77 4- 7j. 
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Figure 2: The individual contributions to the cross section from events where both photons are generated according to 
the PYTHIA diphoton matrix element and from events where one photon originates from initial or final state radiation, 
as functions of the diphoton mass (left), transverse momentum (middle) and azimuthal difference (right). 



Table I: The total diphoton production cross section obtained from the measurement and from the theoretical calculations. 
The PYTHIA 77 calculation involves only the qq — ^ 77 and gg — >■ 77 processes. The pythia 77 + 7] calculation includes 
also the qq — >■ -yyg and gq — >■ 775 processes. 

Cross section (pb) 
Data 12.47 ± 0.21etat ± 3.74,y,t 

RESBOS 11.31±2. 45syst 

DiPHOX 10.58 ± 0.55syst 

PYTHIA 77 + 7j 9.19 
PYTHIA 77 5.03 



4. Results 

The measured total cross section is shown in Table U together with the predictions from the three theoret- 
ical calculations. All three calculations are consistent with the size of the measured cross section within the 
experimental uncertainties. 

Fig. [3] shows the comparison between the measured and predicted diphoton cross sections as functions of 
the diphoton invariant mass M, the diphoton transverse momentum Pt and the difference Acj) between the 
azimuthal angles of the two photons in the event. While the pythia direct calculation (77) fails to describe 
both the scale and shape of the data, including radiation brings the prediction in fair agreement with the data. 
In particular, radiation makes the Pt and At/) distributions harder because of the presence of at least one hard 
jet in the final state of events in which one photon originates from radiation. The mass distributions show a 
reasonable agreement with the data for all predictions above the peak at 30 GeV/c^, particularly in the region 80 
GeV/c^ < M < 150 GeV/c^ relevant to searches for the Higgs boson 1]. However, all predictions underestimate 
the data around and below the peak. In the Pr spectrum all predictions underestimate the data in the region 
between 20 and 50 GeV/c, a feature also observed in the earlier measurements [1, [lB|- For Pt<20 GeV/c, 
where soft gluon resummation is most important, only the resbos prediction describes the data. Discrepancies 
between data and theory are most prominent in the comparison of the measured and predicted distributions of 
A(p. In this case all three predictions fail to describe the data across the whole spectrum. Approaching Aip — tt, 
where soft gluon processes are expected to manifest, the resbos prediction agrees better with the data. In the 
range 1.4 rad<A(^<2.2 rad only the pythia prediction describes the data and remains closest to the data down 
to 1 rad. In the low Acj) tail, which corresponds to the region of low M (<50 GeV/c^), all three predictions are 
lower than the data, although the DIPHOX prediction, by explictly including non-perturbative fragmentation, 
lies closer to the data for A(j)<l rad. 

Fig. 2] shows the comparison between the measured and predicted diphoton cross sections as functions of 
the diphoton rapidity Yy^, the cosine of the polar angle in the Collins-Soper frame cos 9 |27|], and the ratio z 
of the subleading photon Et to leading photon Et- All three predictions agree fairly well with the measured 
da/dY^y, da/dcosd and da/dz, within uncertainties. Exceptions are the predictions of all three calculations 
underestimating the data in the two ends of the cosO spectrum, where again gluon scattering processes and 
associated fragmentation are expected to dominate (T]| . 

In general, all three calculations reproduce most of the main features of the data, as observed in the earlier 
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diphoton cross section measurements |16l . Il7| . However, depending on their approximations, they display 
differences with each other and with the data in certain kinematic regions. There is a problem common to all 
three calculations in the description of events with very low diphoton mass, low azimuthal distance and diphoton 
transverse momentum in the region of the so-called "Guillet shoulder" (20 GeV/c<PT<50 GeV/c) [|. Such 
events include fragmentation at a relatively high rate. The pythia 77 calculation fails completely to describe 
the data both in the scale, where it is low by a factor of 2.5 (see Table |T]), and in the shape, particularly of the 
Ptj ^0 and z distributions, where it predicts a much softer spectrum than the data. This is in agreement with 
the conclusion of Ref. [l^ which tested only pythia 77 as a parton showering MC prediction. 



5. Summary 

In summary, the diphoton production cross section, differential in kinematic variables sensitive to the reaction 
mechanism, is measured using data corresponding to an integrated luminosity of 5.36 fb~^ collected with the 
CDF II detector. The high statistics of the measured sample allows for a higher precision scan over a much more 
extended phase space than previous measurements. The overall systematic uncertainty is limited to about 30%. 
The results of the measurement are compared with three state-of-the-art calculations, applying complementary 
techniques in describing the reaction. All three calculations, within their known limitations, reproduce the main 
features of the data, but none of them describes all aspects of the data. The inclusion of photon radiation in 
the initial and final states significantly improves the pythia parton shower calculation (see the left-hand panels 
in Figures El and 3]) , which is suitable for background simulations in searches for a low mass Higgs boson and 
new phenomena. 
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Figure 3: The measured differential cross sections compared with three theoretical predictions discussed in the text. The 
left windows show the absolute comparisons and the right windows show the fractional deviations of the data from the 
theoretical predictions. Fractional deviations for pythia refer to the 77 + 7j calculation. Note that the vertical axis scales 
differ between fractional deviation plots. The comparisons are made as functions of the diphoton mass (top), transverse 
momentum (middle) and azimuthal angle difference (bottom). The shaded area around the data points indicates the 
total systematic uncertainty of the measurement. 
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Figure 4: The measured differential cross sections compared with three theoretical predictions discussed in the text. The 
left windows show the absolute comparisons and the right windows show the fractional deviations of the data from the 
theoretical predictions. Fractional deviations for pythia refer to the 77 + 7] calculation. Note that the vertical axis 
scales differ between fractional deviation plots. The comparisons are made as functions of the diphoton rapifity (top), 
cosine of the polar angle in the CoUins-Soper frame (middle) and ratio of the subleading photon Et to leading photon 
Et (bottom). The shaded area around the data points indicates the total systematic uncertainty of the measurement. 



